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ABSTRACT 

We present the most energetic BALQSO outflow measured to date, with a 
kinetic luminosity of at least 10^^ ergs s~^, which is 5% of the bolometric lumi- 
nosity of this high Eddington ratio quasar. The associated mass flow rate is 400 
solar masses per year. Such kinetic luminosity and mass flow rate should provide 
strong AGN feedback effects. The outflow is located at about 300 pc from the 
quasar and has a velocity of roughly 8000 km s~^. Our distance and energetic 
measurements are based in large part on the identiflcation and measurement of 
S IV and S IV* BALs. The use of this high ionization species allows us to general- 
ize the result to the majority of high ionization BALQSOs that are identifled by 
their C IV absorption. We also report the energetics of two other outflows seen 
in another object using the same technique. The distances of all 3 outflows from 
the central source (100-2000pc) suggest that we observe BAL troughs much far- 
ther away from the central source than the assumed acceleration region of these 
outflows (0.01-O.lpc). 

Subject headings: galaxies: quasars — galaxies: individual (SDSS J1106-I-1939, SDSS 
J1512-I-1119) — line: formation — quasars: absorption lines 
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1. INTRODUCTION 

Broad absorption line (BAL) outflows are observed as blueshifted troughs in the 



rest- 



2008 



rame spectrum of ^- ^20 % of quasars ( iHewett &: Foltz 



2003 



Knigge et al 



Ganguly fc Brotherton 



20081). The energy, mass, and momentum carried by these outflows 



are t 



lought to play a crucial role in shaping the early universe and 



(e.g. 


Scannapieco & Oh 


2004 


2009; 


Ciotti et al. 


2009, 


2010 



Levine & Gnedin 


2005 


Ostriker et al.I 


|201o|). ' 



universe 


and ( 


ictating its evolution 


Hopkins 


et al. 


20061 


Cattaneo et al. 



deduced from the detection rate of these mass outflows, allows for efficient interaction with 
the surrounding medium, and these outflows carry thousands of times more mass flux per 
unit of kinetic luminosity than the collimated relativistic jets observed in 5 - 10% of all 
AGNs. Theoretical studies and simulations show that this so-called AGN feedback can 
provide an explanation for a variety of observations, from the chemical enrichment of the 



intergalactic medium, to the self regulation of t 



and of the galactic 



Hopkins et al. 



j ulge (e.g. 



2009 



Elvis 



Silk fc Rees 



2006 



1998 



le growth of the 



Di Matteo et al 



supermassive black- 



lole 



2005 



Germain et al. 



2009 



and references therein). 



Quantifying AGN feedback requires estimating the kinetic luminosity [Ek) and 
mass-flow rate (M) of the outflows. These quantities can be computed in cases where we 
are able to estimate the d istance R to the ou tflowing material from the central source 



(see Equations 6 and 7 in 



Borguet et al. 



2OI2I ). Following the deflnition of the ionization 



parameter of the plasma Uu cx I/uhR? where Uh is the hydrogen number density, 
the distance can be obtained for outflows of known ionization state and density (see 
elaboration in Section |5]). The research program developed by our team ha s led to the 



determination of M and Eh in several 



2010 



Bautista et al. 



2010 



Aoki et al. 



quasar outflows (e.g. iMoe et al. 



2011 



Borguet et al. 



2009 



Dunn et al. 



20121 ) by using the population 



^Based on observations collected at the European Southern Observatory, Chile, PID:87.B- 
0229 
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ratio of coUisionaly excited state s to the resonance leve ls of singly ionized species (e.g. Si ll, 



Fell) as density diagnostics (see 



Crenshaw et al. 



2003. for a review). The lower detection 



rate of these low ionization outflows in spectroscopic surveys raises the question of whether 



the determinations obtained for these objects are 
ionization C IV broad absorption line quasars (see 



representative o 



Dunn et al. 



the ubiquitous high 



20121 . hereafter Paper I). 



One way to alleviate this uncertainty is to target objects which possess absorption 
troughs from excited states of high ionization species, where Siv/Siv* AA1062.66, 1072.97 
are especially promising. These transitions appear at wavelengths blueward of Lya and 
therefore suffer from blending with the Lyman forest in high redshift objects. However, 
since the ratio of C iv and S iv ionic fractions as a function of the ionization parameter (Uh) 
is relatively constant (see paper I) they must arise from the same photoionized plasma. 
The ionization similarity of Civ and Siv makes Siv/Siv* outflows a much better agent 
than the usual low ionization species for the determination of the feedback from the high 
ionization outflows (see discussion in Section [5^ . We consequently developed a research 
program that measured the sample properties of these S iv outflows (Paper I) and analyzed 
the photoionization and chemical abundances of one such outflow (Borguet et al. 2012b, 
hereafter paper II). In this paper, we present the analysis of VLT/X-shooter spectra of 
two SDSS BAL quasars, SDSS J1106+1939 and SDSS J1512+1119, for which the presence 
of Siv/Siv* troughs allows us to place constraints on the location and energetics of the 
outflow in SDSS J1106+1939 and of two separate outflows in SDSS J1512+1119. 

The plan of the paper is as follows: In § [2] we present the VLT/X-shooter observations 
of SDSS J1106+1939 and SDSS J1512+1119 along with the reduction of the data. In § [3] 
we identify the spectral features associated with S IV and S iv* and measure the column 
densities for various ionic species associated with each outflow. Photoionization models are 
crucial for finding the total hydrogen column density of the outflows and their ionization 
equilibrium. We present these models in § 4. In § [5] we derive the parameters necessary to 
determine R, Ek and M of the S iv outflows. In § [6l we evaluate the robustness of each step 
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in the analysis and conclude that Ek = 10^^ ergs s ^ is a firm lower limit. We discuss our 
results in § [71 

2. OBSERVATIONS AND DATA REDUCTION 

We selected our targets SDSS J1106+1939 (J2000: RA=11 06 45.05, dec=+19 39 
29.1, g=19.4) and SDSS J1512+1119 (J2000: RA=15 12 49.29, dec=+ll 19 29.4, g=17.7) 
from an extensive flux limited search of the SDSS DR7 catalog for objects indicative of the 
presence of S iv absorption on the basis of their SDSS spectrum. We observed these two 
objects with VLT/X-shooter as part of our program 87.B-0229 (PI: Benn) in April 2011 and 
March 2012, respectively. X-shooter is the second generation, medium spectr al resolution 



(R ^ 6000-9000) spectrograph installed at the Cassegrain focus of VLT/UT2 (IVernet et al. 



20111 ). The unique design of the instrument, in which the incoming light is split into three 
independent arms (UVB, VIS and NIR) each composed of a prism-cross-dispersed echelle 
spectrograph, allows the simultaneous covering of a wide spectral band (3000 A to 24000 
A) in a single exposure. The total integration time for each object in the UVB, VIS and 
NIR arms were 8400, 8400, and 8700 s, respectively. 

We reduced the SDSS J1106+1939 spectra in an identical fashion to the one of 



SDSS J1512+1119 (detailed in Paper II): we rectifie d and wavelength ca. 



ibrated the two 



201l[ ). then extracted 



dimensional spectra using the ESQ Reflex workflow (iBallester et al. 
one dimensional spectra using an optimal extraction algorithm and finally flux calibrated 
the resulting data with the spectroscopic observations of a standard star observed the same 
day as the quasar. The resulting flux calibrated spectrum of SDSS Jl 106+1939 is presented 
in Figure [TJ We only present the UVB+VIS spectrum and present the few additional 
diagnostic absorption lines detected in the NIR portion of the spectrum in Fig. |2j 
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Fig. 1. — Reduced and fluxed UVB+VIS X-shooter spectra of tlie quasar SDSS J1106+1939. 
We indicate tlie position of the absorption lines associated with the intrinsic outflow. The 
VIS part of the spectrum {Xobs > 5700 A) has not been corrected for atmospheric absorption. 
This does not affect our study since our diagnostic lines are located in regions free of such 
contamination. The dashed line represents our unabsorbed emission model (see Section [STT]) . 
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SPECTRAL FITTING 



The radial velocity values across t he absorption troughs are determined with respect 



to the systemic redshift of the quasar. 



Hewett &: WildI (I2OIOI ) report an improved redshift 



valuqjfor these two SDSS DR7 quasars by cross correlating a quasar emission line template 
to the observed SDSS spectrum. These improved redshifts are z = 3.0377 ± 0.0021 for 
SDSS J1106+1939 and z = 2.1062 ± 0.0020 for SDSS J1512+1119. Note that while the NIR 
range of the X-shooter observations of J1106+1939 covers the H/3+ [Olll] A5007 and also 
the Mgll rest frame regions, these portions of the spectrum are located within dense and 
strong H2O and CO2 atmospheric bands preventing us from determining a more accurate 
redshift from the fit of these emission features. We also examined the expected position of 
[0 11] A3727 and found no emission feature in its spectral vicinity. 



3.1. Unabsorbed Emission Model 

Deriving ionic column densities from absorption troughs requires knowledge of the 
underlying unabsorbed emission Fo(A). In AGNs the UV unabsorbed emission source can 
generally be decomposed into two components: a continuum source described by a power 
law and emission lines usually modeled by Gaussian profiles that are divided into broad 
emission lines (BEL, full width at half maximum (FWHM) > 1000 km s~^) and narrow 
emission lines (NEL, FWHM < 1000 km s"^). 

Due to the presence of strong absorption shortward of the rest-frame Lya emission 
in SDSS J1106+1939, we fit the continuum in the range 3200 - 9000 A using a single 
power law of the form F{\) = Fiioo(A/1100)", where Fnoo is the observed flux at 1100 A 
in the rest-frame of the object, on regions longward of the Lya emission suspected free 
of absorption or emission lines. After correcting the spectrum for the galactic extinction 



^Available online at: 



http: / /das. sdss.org/va/Hewett_Wild_dr7qso_newz/ 
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(E(B-V) = 0.025. ISchlegel et al.l 119981 ) using the reddening curve of ICardelli et al.l (119891 ) 
we find -Fnoo = 12.8 x 10^^^ ergs/s/A/cm^ and a ~ —0.95. We fit the prominent Lya and 
C IV emission hues using a sum of Gaussians and a sphne fit for the Si lv+0 iv] blend. We 
fit the Fell and other weak emission hues longward of 1600 A in the rest frame using a 
phenomenological sphne fit. In the far UV, we build an O vi emission model by scaling and 
shifting to the proper position the C iv emission line model. The scaling factor is chosen 
to be the smallest such as no emission remains above the constructed emission model. 
This approach could strongly underestimate the true underlying O VI emission. We fit the 
NIR unabsorbed emission using a smooth spline fit. The resulting emission model for the 
UVB+VIS region is shown in Figure [H Note that the single power law model for the whole 
UVB+VIS range most likely provides an overestimation of the true underlying continuum 
shortward of the O yi emission, a region in which a softe r power law is usu ally adopted 



(IKorista et al. 



1992 



Zheng et al 



1997 



Arav et al.n2001b 



Teller et al 



20021 ). However, 



modeling of the unabsorbed emission in that region is not important for our BAL analysis 
due to the fact that it contains only heavily blended diagnostic lines coupled with a limited 
signal to noise ratio (S/N). 

Given the overall high S/N of the SDSS J1512+1119 data spectrum (S/N - SO- 
TO over most of the UVB/VIS range) and absence of wide absorption troughs, we fit the 
unabsorbed emission lines us ing a smooth third o rder spline fit after fitting the derredenned 

Schlegel et al.lll998[ ) with a power law with Fnoo = 77.0 x 10~^^ 



continuum (E(B-V) = 0.051, 
ergs/s/A/cm^ and a ~ —0.403 (see Paper II for details). 



3.2. BAL Column Density Measurements of SDSS J1106+1939 



The X-shooter spectrum of SDSS J1106+1939 exhibits wide absorption troughs 
associated with Hi, CiV, Siiv, Nv, Ovi, Siv/Siv* and Pv ionic species. The CiV 



absorption trough (and most oth er troughs as we 



broad absorption line (BAL, see 



Weymann et al. 



I) sat isfies the observational definition of a 



19911 ). In addition, we identify absorption 
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at the same kinematic location associated with He I* and the lower ionization species Mg ll, 
C II, Al II, and Al ill (see Figures [U and [2]) . 

The column density of an ionic species i associated with a given kinematic component 
is estimated by modeling the residual intensity /j(A) = FobsW / Fo{\) as a function of the 
radial velocity. The simplest modeling technique is the apparent optical depth (AOD) 
method where Ti{v) = — ln(/j(f)), which is then converted to a column density Ni(v) us ing 



the appropriate aton iic/physica. 1 constants (see Eq u ation 9 in 



However, our group (lArav 



Gabel et al 



Churchill et al. 



19971 : 



2005b) and others ( 



1999 



Arav et al.i 



Bar 



Ganguly et al 



ow et al 



999al. 



1997 



b. 1200 la 



Savage fc Sembach 



2002 



Hamann et al 



2003 



1997 



1991) 



Scott et a. 



Telfer et al. 



2004 



1998 



19991 ) showed that column densities derived from the 
apparent optical depth analysis of BAL troughs are unreliable due to non-black saturation 
in the troughs. In particular, in Paper II we showed that the true C iv optical depth in the 
deepest outflow component of SDSS J1512+1119, is ~1000 times greater in the core of the 
absorption profile than the value deduced from the AOD method. 

To account for non-black saturation in unblended doublets or multiplet troughs 



from the same ion, we routinely use the partia 



absorption models (e.g. 



Arav et al. 



1999a, 



2008 



covering (PC 



Edmonds et al. 



and power-law (PL 



2011 



Borguet et al. 



2012 . 



for details). However, the intrinsic width of most absorption troughs in the spectrum of 
SDSS J1106+1939 causes self-blending of troughs from most of the observed doublets (CiV, 
Nv, Ovi, Mgll, Aim, Siiv and Pv). As a result, the pure PC and PL methods cannot 



the s t udy of BAL qua sar spectra (e.g. 



Korista et al. 


1992; 


Arav et al. 


1999b: 


de Kool et al. 



2002 



Moe et al. 



2009 



and references therein). The main assumption made when using 



this technique is that the physical properties of the absorbing gas do not signif icantly 



change as a function of the radial velocity for a given kinematic component (e.g. 



2009 



Dunn et al. 



Moe et al. 



20101 ). Therefore, the optical depth profile as a function of velocity 



is assumed to be proportional for all lines so that a single scaled template is used to 
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r epro d uce the observed absorption troughs. 



tested (IKorista et aL 



2008 



Moe et aL 



2009 



n the three cases where this assumption was 



Dunn et aL 



20101), it was shown to be well 



consistent with the data. As we will show here, we obtain good fits to most troughs in SDSS 
J1106+1939 using this simple and restrictive assumption. During the fitting procedure, we 
assume that doublet lines are not affected by non-black saturation, as a first approximation, 
and relax that constraint when needed in order to obtain a better fit to the observed line 
profile. It is important to note that in most doublet and multiplet cases, template fitting 
can indicate whether the trough is highly s aturated or whet her the actual column density 
is close to the non-saturated AOD case (see 



Moe et al. 



2001 for details). 



We choose the P v AA 1117.98,1128.01 doublet lines as a starting point for the template 
building. Our main motivations in the choice of P V are the high velocity separation between 
its components {Av ~ 2700 km s^^), and the absence of other strong intrinsic absorption 
and emission lines in that region of the spectrum. The first step in our template building 
process consists of matching the blue wing profile of the P v BAL since that part of the 
trough is not blended with kinematic components from the red line of the doublet. The 
absorption structure visible in the P V blue wing suggests that more than one kinematic 
component might be present for the blue transition. Therefore, for our first model we choose 
the optical depth template of each line to be a sum of two Gaussians Gl and G2 (Model 1). 

In Figure |2] we show the result of the best fit to the moderately blended transitions 
in the SDSS J1106-I-1939 spectrum using Model 1. The FWHM [wi) and central positions 
(f j) of the two Gaussians have been fit simultaneously to reproduce the P v BAL troughs. 
The fit suggests that the blue wing of the P v BAL profile is well represented by a simple 
model in which the first Gaussian Gl is characterized by wi = 1150 km s^^, vi = —8250 km 
s~^ and the second Gaussian G2 by W2 = 2400 km s~^, V2 = —10100 km s^^. An absolute 
lower limit on the P v column density present in Gl and G2 is placed by using the AOD 
technique on the blue transition whose profile is known and unblended. Relaxing the 2:1 
constraint of optical depth between the doublet lines of P V, we are able to obtain a better 
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match of the observed hne profile at the location of the core of Gl PvA1128.01 (see first 
panel of Figure [2]) by increasing the optical depth associated with that transition. We place 
an upper limit on P V column density in Gl by using the PC technique on this 1:1.1 optical 
depth ratio doublet (see Tabled]). This estimation constitutes a hard upper limit since 
the red P v line of that component is likely blended by the P V blue transition of a lower 
velocity system (see Model 2) such that the true ratio of optical depth for that system is 
within 1:2 to 1:1.1. 

The optical depth model derived from the P V BAL fit is also able to reproduce a 
significant part of the blended Siv/Siv* BAL. The Siv/Siv* model presented in Figure [2] 
considers the presence of the resonance A1062.97 and excited A1073.01 line (the oscillator 
strength weighted mean of the blend of the two close-by Eiow = 951 cm~^ S IV transitions, 
see paper II for details) for which Av ~ 2900 km s~^. Due to the uncertainty in the 
Ovi emission contribution in the blue wing of S iv (see Section [3TT]) . the column density 
associated with component G2 for that transition is less certain and could be underestimated 
by up to a factor of two if we would not have scaled down the shifted C iv emission model. 
Model 1 is also able to reproduce the absorption profile observed in the other species such 
as He I*, Mgll, Cll, Alll, and Aim as can be seen in the remaining panels of Figure [2] in 
which the maximum optical depths are the only free parameters of the models. The fit 
of these profiles was performed without having to relax the 2:1 assumption for the ratio 
of optical depth between the strongest and weakest doublet components, consistent with 
non saturated troughs. We report the column densities derived for each ionic species in 
Table [H In He I*, we do not detect any absorption related to component G2 and report an 
upper limit on the column density for that component by scaling the Gaussian template 
assuming that that the noise level could hide a 2a detection for the strongest He I* A3889.80 
transition. The Civ, Siiv, Ovi, Nv doublet are both more heavily self-blended than 
P V and S IV because of the smaller velocity separation between the doublet components, 
and are much more saturated due to their higher elemental abundance. This prevents 
us from obtaining reliable fits to their absorption troughs, and only lower limits on their 
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column densities can be derived. Blending of the Nv, Hi and Sim troughs prevents us 
from deriving useful column densities for these species as well (beside the fact that N v, H I 
must be heavily saturated). We however note that the lower limits derived from the AOD 
technique for these species are consistent with the predicted values in our photoionization 
models. We place an upper limit on the column density of the non detected Fe ll by scaling 
the template of each Gaussian to within 2a at the expected location of the Fell A1608.45 
troughs. We did not use the stronger Fell A2600.17 that is located in a region of poor SN 
at the blue edge of the NIR detector or the Fell A2382.76 that is severely blended by a H2O 
atmospheric band. 

While Model 1 provides a reasonable fit to the observed absorption profile associated 
with most of the ionic species, it fails to reproduce the red wing observed in the BAL 
profile of high ionization species such as Pv, Siv+ Siv* or Siiv (see Figure [2]). This 
observation suggests the existence of an additional kinematic component associated with 
the high ionization species. In order to better fit the red wing of the high ionization lines, 
we built a second template. Model 2. This second optical depth template is also composed 
of two functions: a high velocity Gaussian F2 with identical parameters as G2, and Fl: 
a low velocity modified Gaussian profile. The result of the simultaneous fit, performed 
with the same constraints as with Model 1, is presented for Pv, SlV+ Siv* and Si IV in 
Figure [3], and the associated column densities are shown in Table [21 Here the overall fit 
of the P V BAL is much better, as well as the fit to the S iv+S iv* profile. The reported 
column densities are consistent within ~30% with the one derived using Model 1. While 
Model 2 provides a better fit to the high ionization P v and S iv+ S iv* BAL (but also 
Silv), the lack of constraints on the shape of the low velocity component template makes 
the derived column densities too model dependent to be reliable. Moreover, the small 
fractional difference in column density from that component does not significantly affect 
the photoionization modeling since the bulk of the computed column density is located in 
the two main components Gl and G2. For these reasons, we will use the column densities 
reported in Table [H 
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Fig. 2. — Fits to the mildly blended high- ionization absorption troughs (P v, S iv and Si iv), 
as well as the non saturated troughs (He I*, Mg ll, C ll, Al ll, Al ill) in the X-shooter spectrum 
of SDSS J1106+1939. Each transition is fitted with an identical optical depth profile whose 
shape is described by a sum of two Gaussians G1+G2 (Model 1, see Section I3.2p leaving 
only the maximum optical depths as a free parameter in each fit. Gl is represented by a 
dotted line and G2 by a dashed line. The full model for the trough is plotted as a solid line. 
This simple model provides a good fit to He I* and the low ionization species as well as to 
the unblended blue wing of the P v absorption. The poor fit of the P v, S iv+ S iv* and 
Si IV BAL red wings suggests the existence of another component that is not seen in the low 
ionization species (see text and Figure [3]). 
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Table 1. Ionic column densities in the outflow of SDSS J1106+1939 using Model 1 



Ion 


Gl (-8250 km s^i) 


G2 (-10100 km s^i) 


A r\ r-\T~\i-a,r\ ^ 




AOD 


PC 


AOD 


PC 


(CA+G2) 






(10^2 cm-2) 


(10^2 cm-2) 


(10^2 cm-2) 


('1012 f.m-21 


He I* 


450lfo 


48012™ 


< 140 


... 


480l3« 


Cii 






2201^° 




730lf„ 


Cii* 


9801^20 




5301^° 




15001^° 


Alii 


161} 




< 7.1 




161? 


Aim 


2401}° 


2501}° 


1801}° 


1801}^ 


4301?° 


Pv 


1300130 


45001?™ 


19001^° 


27001^™° 


3200-9900^ 


Siv 


380001^2™° - 




52001}°° 




440001??°™ 


Siv* 


leoool^oo 




25001}°° 




190001™° 


Mgii 


14011° 


1401}° 


841^5 


85lf 


2901?° 


Fell 


< 160 




< 370 




< 540 



a) Column density values adopted for photofonizatfon modeling (see Section HJ. 
b) The range translates the uncertainty on the P v column density due to the possible non black saturation 

affecting component Gl (see text), 
c) The upper error for the black Siv Gl component has been computed by finding the maximum optical 
depth that can be hidden until the wings of component Gl become too broad for the observed line profile. 



Table 2. Ionic column densities in the outflow of SDSS J1106-f 1939 using Model 2 



fon Fl (-8250 km s'^) 

AOD PC 
(10^2 cm-2) (^0^2 cm-2) 



F2 (-10100 km s-^) 
AOD PC 

(10^2 cm-2) (^Ql2 ^^^-2) 



Pv 



3300 



+40 
-30 



5700 



+3200 
-110 



-1+4200 



S IV 470001^200 



Siv* 



2600013™ 



1800 
5500 



+30 
-30 

-560 
-140 



2400; 



-320 
-30 



29001? 



+530 
30 
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We presented two phenomenological decomposition of the BAL profiles that allows 
us to estimate the column densities associated with S iv and S iv* to better than within 
50% (see Tabled] and [2]). The larger uncertainty affecting the Siv column density in the 
resonance state has been estimated by scaling the Gl template so that the wings of the 
model are no longer consistent with the observed profile. In the last column of Table [Tj 
we report the adopted column densities for the photoionization modeling, as well as the 
statistical error affecting the measurements (taking only the photon noise into account). 
For photoionization modeling we choose to sum components Gl and G2 for the following 
reason: The ratio of almost all ionic column density of the same ions is roughly 2:1 between 
Gl and G2, suggesting very similar photoionization solution. The exception is Siv where 
the ratio is ~ 7. However as we noted above, this is mainly due to our choice of a minimal 
O VI BEL. If we choose an O vi BEL as strong as the C iv BEL the column density of S iv 
doubles and triples if we choose a 1.5 times stronger O vi BEL. This will make the S iv ratio 
between components Gl and G2 ~3.5 and ~ 2 respectively. We therefore conclude that 
it is highly plausible that both components arise from the same outflow and have similar 
ionization equilibria, which justify adding them together for the purpose of photoionization 
modelling and the extraction of the physical parameters of the outflow. When available, 
we choose to use the value reported in the PC column as the measurement. If only AOD 
determination is available, we will consider the reported value minus the error as a lower 
limit during the photoionization analysis since conceptually no information about non-black 
saturation effects can be obtained from singlet lines. For Pv, due to uncertainty in the 
column density associated with component Gl, we report the range of values that are 
allowed within the absolute lower and upper limits placed on the column density. 

The template decomposition also demonstrates the detection of a S iv* BAL outflow 
and the fact that the column density associated with the excited state is lower than the one 
measured in the resonance level. The detection of the S iv* BAL is secured by the good 
kinematic match of the template fit. Furthermore, if that part of the profile would have 
been related to another low velocity S iv system rather than S iv*, then no kinematic match 
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would be observed with not only P v, but also with higher abundance species like C IV, O vi 
or Nv that would have likely been heavily saturated. 

We note that the observation of absorption associated with the low ionization Al ill and 
Sim lines suggests the presence of absorption due to the FelllA 1122.524 blending with the 
P V BAL line. Such a blend could affect the P v column density derived during our template 
fitting procedure if the Fe ill column density was significant enough. The photoionization 
model presented for the SDSS J1106+1939 outflow in Section H] predict an optical depth 
r < 0.1 for that Felll line, which does not affect our P v column density determination. 

3.3. Column Density in SDSS J1512+1119 

Our second target, SDSS J1512+1119, displays a Civ absorption trough whose width 
is close to a BAL {Av ~ 1700 km s~^, see Paper II). However, the Si IV lines are narrow 
enough and unblended {Av ~ 2000 km s~^) to be adopted as a template to identify the 
kinematic components in other ionic species. Using that template, we reported in Paper II 
the existence of five main kinematic components with centroids at velocities of f ~ -2100, 
-1850, -1500, -1050, and -520 km s~^ (labeled component 1-5). In this paper we concentrate 
on determining the distance and energetics associated with the two main components: 
component 2 and 4. As shown in Figure H] in which we use the Siiv A1393.76 line as 
template, absorption troughs associated with S IV and S iv* are observed in kinematic 
component 2 (f ~ -1850 ^^) and Siv absorption is also detected in kinematic component 
4 (t; ~ -1050 ^^). The in depth analysis of the SiV and SiV* troughs performed in Paper 
II revealed that the blend affecting the red wing of the component 2 in S iv* A1072.97 is 
identified as the ten times weaker transition of SiV* A1073.51. This allowed us to derive 
rehable estimates of N{S iv) ~ WGOOttmo x 10^^ cm'^ and N{S IV*) ~ 16000tf^|||] x lO^^ 
cm~^ for that component. 

For component 4, we estimate the S IV column density by using the AOD model on 



-17- 




1070 1080 1090 1100 
Rest Wavelength (A) 



1110 




1020 1030 1040 1050 1060 
Rest Wavelength (A) 




1320 1330 1340 1350 1360 1370 1380 1390 
Rest Wavelength (A) 



Fig. 3. — Fitting the high ionization Pv, Siv+Siv* and Siiv BAL profiles with the Model 
2 optical depth profile. As for Model 1, the optical depth distribution associated with 
a transition is assumed to be known, leaving only the maximum optical depths of each 
components as a free parameter in each fit. Model 2 is composed of a sum of two functions: 
Fl (a Gaussian with a modified red wing) represented by a dotted line and F2 (a Gaussian 
with identical parameters as G2) represented by a dashed line. The total model for the 
trough is plotted as a solid line. For P v the fit requires a ratio of optical depth between 
both components of the doublets of 1:1.1, suggesting the possibility of moderate non-black 
saturation. 
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the normalized line profile (i.e. t{v) = — ln(/j(f))), converting it to column densities and 
integrate it over the trough to obtain a lower limit of A^(Slv) > 263 x 10^^ cm~^. We do 
not detect an absorption trough associated with Siv* for that component, so we report an 
upper limit on the column density for that transition by scaling the S iv template, assuming 
that the noise could hide up to a 2 a detection and find A^(Slv*) < 25 x 10^^ cm^^. Due 
to the narrowness of the line profile and absence of significant blending, we are able to 
measure the Si iv column density and place limits on the ionic column for other species that 
we report in Table [31 The constraints derived from the C iv and O VI column densities are 
consistent with the one obtained from the Nv and are not reported in the table. We finally 
place an upper limit on the Pv column density of A^(P v)< 8.1 x 10^^ cm^^ due to its non 
detection in that component. 



4. PHOTOIONIZATION ANALYSIS 



We use photoionization models in order to determine the ionization equilibrium of 
the outflow, its total hydrogen column density (Nh), and to constrain its metallicity. The 
ionization parameter 



H 



Qh 



(where Qh is the source emission rate of hydrogen ionizing photons, R is the distance 
to the absorber from the source, c is the speed of light, and is the hydrogen number 
density) and A^^^ of the outflow are determined by self-consistently solving the ionization 
and thermal bal ance equations with version c08.00 of the spectral synthesis code Cloudy, 



1^ 



last described in 



Ferland et al 



( Il998l ). We assume a plane-p arallel geometry for a gas of 



Lodders et al 



fl2009[ ). Given 



constant ra^ and initially choose solar abundances as given in 
the lack of observational constraints in wavebands outside the X-shooter spectral range for 
both objects, we choose the UV-soft spect ral energy d i stribu tion (SED) model for high 



luminosity radio quiet quasars described in 



Dunn et al 



mm . The use of this SED mod el 



lacking the so-called big blue bump from the classical MF87 SED (IMathews fc Ferland 
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Fig. 4. — Identification of tlie S IV and S IV* intrinsic absorption trouglis in tlie spectrum 
of SDSS J1512+1119. The red wing of tlie trouglis associated with component 2 of Siv* 
A1072.97 is mainly affected by a blend due to the S iv* A1073.51 transition of component 2, 
that allowed us to pinpoint the SiV* column density for that component (see Paper II). 
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Telfer et al 



Dunn et al 



fl2002h over 
for a 



2010 



19871 ). is motivated by the rather soft FUV slopes observed by 
a large sample of HST spectra of typical radio quiet quasars (see 
detailed discussion). Using this SED, we generate a grid of models by varying Nh and 
Uh- Ionic column densities predicted by the models are tabulated and compared with the 



measured values in order to determine the models that best reproduce the data. 



4.1. SDSS J1106+1939 

We compare the ionic column densities given in the last column of Table [1] (see discussion 
in Section [3l2|) with predictions of photoionization models, and use = 10'^ cm~^, which 
we measure for the outflow using the ratio of Siv/Siv* ionic column densities (see 
Section fS.l.ip . We focus our modeling on the troughs from the high ionization species He I*, 
P V, and S IV, as these are the dominant species in the outflow, and S iv is the species that 
allows us to constrain R (see Section lB.l.ll) . We treat He I* as a high ionizatio n species sin ce 



its con centration depends linearly on the fraction of Hell in the gas (see e.g. 



Arav et al. 



2001ah . 



In Figure [5], contours where model predictions match the measured ionic column 
densities are plotted in the Nh — Uh plane. Note that we use the average of the range of 
values given for P v in component Gl and include the upper and lower limits in the error to 
reflect the uncertainty on the exact column density in that component (see Section [X^ . The 
best-fit model is parametrized by logiUn) = 0.0 ± 0.2, and logiNn) = 22.8 ± 0.2 cm'^. The 
errors are determined by the region where P V and S IV bands cross and are correlated, with 
higher ionization parameters corresponding to higher total column densities. Comparison 
of measured and predicted ionic column densities are given in Table |H As can be seen 
from Table HI the best model predicts the column densities of the high ionization species 
to within a factor of 2. The poor fit to ClI and Aim is not physically troubling for the 
following reasons: First, as stated above, we purposefully attempted to find the best fit for 
the dominant high ionization species in the outflows. Second, within the reported error 
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bars for Nh — Uh we can find solutions that yield a much better fit for these two ions. For 
example, a model with \og{UH) = —0.2 and \og{NH) = 22.85 cm~^ produces both the Cll 
and Aim column densities to better than a factor of two, with only a moderate worsening 
of the high ionization lines fit. The drastic changes in the column densities of the singly 
ionized species is d ue to the proximity of the solution to the hydrogen ionization front (see 



Korista et al. 



20081). 



The analysis above relies on the assumption of solar abundances. However, AGN 
outfiows a re known to have moderate supersolar metallicities (e.g. , QSO J2233-606: 



4Z 



0) 



Gabel et al. 



2006 



Mrk 279: Z ~ 2Z, 



0? 



AravetaL 



20071 : SDSS J1512+1119 



Zq < Z < 4:Zq, Paper II). We therefore investigate how sensitive our results are to higher 
metallicity models. Comparing the S IV line on the grid models to that of He I*, we find 
that the abundance of sulphur must be ^ 4 times its solar value. This is because the 
abundance of helium is relatively insensitive to changes in the metallicity of the plasma 
(within ^ 8 — 13% of hydrogen for Z ^ 4:Zq), and He I* is underpredicted compared with 
S IV for higher sulphur abundances. We therefore run a grid of models using a metallicity of 
4 times solar, which is consistent with the results of the works cited above for other AGN 
outflows. Our elemental abundances for t he supersolar metalli city model are determined by 



scaling C, N, O, Mg, Si, Ca, and Fe as in 



Ballero et al 



( 120081 ) while the remaining metals 



are scaled as in Cloudy starburst models. The computed models are shown in Figure |6] 
in which the best-fit model for Z = AZq is characterized by log(f//f) = — 0.5to2; and 
log(A'^i^) = 22.ltai cm~^. The quoted errors on Uh and Nh are determined in a similar 
fashion to the solar metallicity case. We consider this model to be the most physically 
plausible for SDSS J1106+1939, while at the same time it provides a more conservative 
(lower) estimate for the kinetic luminosity of this outfiow (see Section [5]). In Section W?2\ we 
explore the sensitivity of the photoionization solution to a different SED. 
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Table 3. Ionic column densities associated with component 4 
of the SDSS J1512+1119 outflow. 



Ion AOD PC 



(1012 cm-2) (1012 cm-2) 



Hi 


> 1427"* 




Nv 


> 1270 




Sim 


> 5.7 




Si IV 






Siv 


> 263 




Siv* 


< 25 





a) Obtained through AOD modeling of the Ly/3 absorption profile. 



Table 4. Results of the photoionization modeling for SDSS J1106+1939 



Model Zq 4 Zq 

Ion log (^) log (^) 



He I* 

Cii 

Aln 

Aim 

Mgii 

Pv 

Siv 



0.37 
-1.83 
-1.81 
-0.58 
-0.09 
-0.04 
-0.03 



0.08 
-1.32 
-0.74 
-0.02 
0.11 
-0.09 
0.34 
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Fig. 5. — Photoionization modeling of the outflow in SDSS J1106+1939 assuming solar 
abundances. The shaded regions represents the locus of points {Uh,Nh) able to reproduce 
the observed ionic column density for a given species (see Table [1]); dashed lines represent an 
upper limit on the column density of the ion; and dotted dashed lines represent a lower limit 
on the column density. The best model is marked by a diamond and the systemic error on 
the solution is represented by the cross. For clarity we do not represent the Al ll line since 
it is co-spatial with the C ll line. 
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Fig. 6. — Photoionization modeling of the outflow in SDSS J1106+1939 assuming a metal- 
licity ZjZQ = 4. The description of the figure is identical to Figure [51 
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4.2. SDSS J1512+1119 

Component 2 of the outflow of SDSS J1512+1119 was analyzed in Paper II. We found 
hgiUn) ^ -0.9+°;} and hgiNn) ^ 21.9+°;} for the UV-soft SED. Using the SED developed 
in Mathews & Ferland (1987), we found the ionization parameter and column density 
dropped by ^ 0.2 and 0.3 dex, respectively. We determined the metallicity of the gas was 
approximately solar with an upper limit of about 4 times solar. 

For component 4 of this outflow, Si IV is the only ion for which we have a column 
density measurement. However, as seen in Figure [TJ the upper and lower limits of other ions 
help constrain the models. Lower limits on N v and S iv along with the upper limit on P v 
constrain the solution to be located within the region interior to these lines in Figure [71 The 
solutions then lie along the the Siiv band within that region, with —1.7 ^ log(f///) ^ — 1.4 
and 19.8 < logiNn) ^ 20.4 (cm'^). 



5. ESTIMATING THE DISTANCE AND ENERGETICS OF THE S iv 

OUTFLOWS 

Assuming that the absorbing material can be described as a thin {AR/R <C 1), 
partially-filled shell, the m ass flow rate (M) an d kinetic luminosity {Ek) of the outflow are 



given by (see discussion in 



Borguet et al 



20121 ): 



M = ATrRnnmpNnv (2) 
Ek = 27TRQ^impNHV^, (3) 

where R is the distance of the outflow from the central source, Q is the global covering 
fraction of the outflow, /i = 1.4 is the mean atomic mass per proton, rup is the mass of 
the proton, Nh is the total hydrogen column density of the absorber, and v is the radial 
velocity of the kinematic component, which is directly derived from the trough's profile. In 
the rest of this section, we detail the determinations of (or constraints on) R and Q that are 
needed for calculating M and E^. 
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Fig. 7. — Photoionization modeling of component 4 of the outflow in SDSS J1512+1119 
assuming solar abundances. As in Figure El a measured column density is represented by a 
shaded area (showing the uncertainties), an upper limit on a column density is represented 
by a dashed line and a lower limit by a dotted dashed line. The region of the phase space 
in which models are able to reproduce the observed constraints is roughly delimited by the 
trapezium. 
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5.1. Determining R: the Distance of the Outflow from the Central Source 

Measuring R is crucial for estimating of M and Ek, and is also essential for 
understanding the relationship of the outflows with the host galaxy and its surroundings. 
In Section HI we derived the ionization parameter (Uh) for each outflow. Therefore, a 
knowledge of Uh and Qh allows us to solve for R directly from Equation ([1]). 



5.1.1. Determining 

In highly ionized plasma the hydrogen number density is related to the electron number 
density through Ug ^ 1.2 n^. Under the assumption of collisional excitation, the ratio of 
level populations between S iv* (E=951 cm^^) and S lv(E=0 cm^^) provides a direct probe 
of He. Using the column densities for Siv and Siv* reported in Table [Hand the electron 
temperature found for our best-flt Cloudy models (10,000 K, a weighted average for SiV 
across the slab), we flnd log(ne) = 4.11;5^q;37 cm~^ for the outflow in SDSS J1106+1939 (see 
Figure IS]). The negative error quoted on the electron number density is computed from the 
uncertainties on the S iv and S IV* column densities given in Table [H In this case the ~ 50% 
statistical error is dominant compared with possible systematic errors due to our continuum 
placement, as well as the assumption of a Gaussian distribution for the absorbing material 
in the trough. However, the positive statistical error on the S IV* to S iv ratio is dominated 
by the 2.5% quoted negative error on the S iv column density (see Table [1]). In this case the 
systematic errors mentioned above should dominate. In order to derive a more physically 
plausible and conservative positive statistical error on the S iv* to S iv ratio we use the 
results of our FH-F2 model (Table |2]). The diflerent distribution of the absorbing material 
for the FH-F2 model causes a much larger deviation in the derived Siv* and Siv column 
density than possible errors due to our continuum placement. We therefore use the higher 
S IV* to S IV ratio of the F1+F2 model as the positive la error for n^. 

Using the Siv and SiV* column densities reported in Table |3] we flnd log(r;,e) < 3.34 
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cm~^ for component 4 of the SDSS J1512+1119 outflow (see Figure [8]). In Paper II we 
showed that no information about n^. could be directly derived from the line profile of S IV 
and Siv* for component 2 of the SDSS J1512+1119 outflow. We therefore used the ratio 
of population between excited states of Cm* A1175 to constrain the electronic density and 
found log(ne) = 5.4^o!6 cm~^. 



5.1.2. Determining Qh and LboI 



We compute Qh as well as the bolometric luminosity LboI for each object by fitting 
the UV-soft SED to the measured flux (corrected for G alactic reddening) at 1100 A (in 



Dunn et al 



f l2010h . We obtain 



the rest-frame) and following the procedure outlined in 
\ogQH = 57.1 s-i and log LboI = 47.2 erg s^^ for SDSS J1106+1939 and logg^ = 57.4 s^^ 
and hglBoi = 47.6 erg s^^ for SDSS J1512+1119. 



5.2. Constraining Vl 



For a full discussion about const raining Q that is n eeded for Equations (2) and (3) 



we refer the reader to Section 5.2 in 



Dunn et al 



([20101) • Here we reproduce the main 



arguments and then concentrate on the case of S iv outflows. There is no direct way to 
obtain the of a given outflow from its spectrum as we only see the material located along 
the line of sight. Therefore, the common procedure is to use the fraction of quasars that 
show outflows as a pr oxy for Q. Statistically, C iv BALs are seen in 20% of all quasars (e.g. 



Hewett fc Folt2 



20031 ). and therefore an Q = 0.2 is used for high ionization outflows. Since 



the ratio of C iv and S iv ionic fraction as a function of Uh is relatively constant (see Paper 
I) they arise from the same photoionized plasma. Therefore, we can simply choose Q = 0.2 
for the S IV outflows. However, to be more conservative we will multiply the Cw Q hj the 
fraction of C iv BALQSOs that show corresponding S iv troughs. The rationale for this is 
that, being a less abundant ion than CiV, S iv may only arise in high outflows and 
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Fig. 8. — Density diagnostic (and errors) for SDSS Jl 106+1939 and kinematic component 
4 of the SDSS J1512+1119 outflow. We plot the theoretical ratio of the level population of 
the first excited states of Sw {E = 951 cm~^) to the level population of the ground state 
versus the electron number density Ug for three representative temperatures. The density 
diagnostic provided is fairly temperature independent for the range of densities represented 
in this plot. We determine the electron number density associated with the two outflows 
using the curve for the temperature of 10,000 K found for our best-fit Cloudy models (see 
text) . 
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therefore a correction factor to that effect may be needed. 



We only compare BAL outflows for two reasons. First, the S iv troughs are located 
within the Lyman forest. The spectral coverage of the Sloan Digital Sky Survey (SDSS) 
can show these Siv/Siv* troughs only for redshifts z > 2.8, where the forest is very 
thick and greatly complicates the identification of narrow S iv troughs. Second, the C IV 
outflows in both SDSS J1106+1939 and SDSS J1512+1119 adhere to the BAL definition 



(IWeymann et al 



199ll ). and therefore the comparison should be made with BALs only. 



We note that nar rower C IV absorption is muc h more prevalent in quasars, up to 60% 



Ganguly &: Brotherton 



20081 ). and therefore their Q should be larger 



of all objects (see 
accordingly. 

The sample we showed in Paper I (their Table 1) contains 24 objects that show 
apprecia ble CiV absorption. F rom these 24 objects, 15 show a Civ BAL trough (again 



using the 



Weymann et al 



1991 



definition). In Paper I, we identified clear cases of S iv/S IV* 
outflows in three of these objects. This was done using the very restrictive criterion of 
matching a Si IV absorption template to both S iv and S iv* troughs (see Paper I, Figure 
1). However, for our purposes here, we need to account for both pure Siv outflow (since 
low density outflows may not show S iv* even when S iv absorption is unambiguous) and 
also cases where the S iv and S IV* are so wide that they blend into one trough. We 
find 3 such strong cases: SDSS J0844+0503 that shows roughly 4000 km s~^ continuous 
Siv/Siv* trough at the expected velocity; SDSS J1051+1532 that shows ~ 6000 km s"^ 
continuous Siv/Siv* trough at the expected velocity; and SDSS J1503+3641 that shows 
a good kinematic match for both S IV and S IV* based on the Si IV absorption template. 
There are also 3 plausible cases of matching narrower Siv and/or Siv* absorption features, 
but we do not include them as the possibility of a false positive identification of unrelated 
Lya absorption features is quite significant. We note that the objects we analyse in this 
paper are not part of the the sample presented in paper I; the SDSS spectrum of SDSS 
J1512+1119 does not cover the Siv/SlV* troughs region, and SDSS J1106+1939 is below 
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the magnitude cut-off of that survey. 

We thus have a total of 6 secure detections of S iv outflow troughs within the 15 C IV 
BAL sample of paper 1 (the three objects shown in Paper I Figure 1, and the three objects 
discussed in the paragraph above). Therefore, a conservative lower limit on the fraction 
of Civ BALQSOs that also show Siv outflows, is 6/15=40%. We will use this correction 
factor for the canonical C iv BALQSOs Q = 0.2, which yields Q = 0.08 for the S iv outflows. 



5.3. Results 

Inserting the parameters derived in Sections H] and Elinto Equations (1), (2) and (3), 
we calculate the distance and energetics for the three S IV outflows. These quantities, as 
well as most of the parameters they are derived from, are listed in Table |5l For outflows 
where ranges of parameters were determined, we show the corresponding ranges in R, M 
and Ek- 

In Table [5] we present 3 different models for the outflow of SDSS J1106+1939: 

• Using solar abundances and the UV-soft SED (see Section H] for details). The 
advantage of this model is its simplicity and the ease of comparison with any other 
outflows that are modelled with solar abundances. 

Using solar abundances an d the canonical MF87 SED developed for AGN by 



Mathews &: Ferlandl ( 119871 ) (see Section [672] below for details), to give a concrete 



example of the results for a different SED. 

• A model with metallicity 4 times higher than solar, combined with the UV-soft SED 
(see Section Wl] for details). As noted in Section Wl] we consider this model to be the 
most physically plausible for SDSS J1106+1939. 



The UV-soft SED we use is the most appropriate given the spectral data on hand, and the 
chemical abundances for Z = AZq give a slightly better fit to the measured column densities 
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than the pure solar case. At the same time, it provides a more conservative (lower) estimate 
for the kinetic luminosity of this outflow. We therefore use this model as the representative 
result for the outflow. 

For the two separate outflow components of SDSS J1512+1119 (C2 and C4), we show 
only the result of models using solar abundances and the UV-soft SED. For comparison, 
we also show in Ta ble [5] the parana eters derived for our previous Ek record holder: 



SDSS J0838+2955 fMoe et al 



2009 



but see the correction by a 0.5 factor reported in 
Edmonds et aLiGOlllV We note that M and Ek for SDSS J0838+2955 were derived using 
Q = 0.2. Therefore, for a fair comparison with the Siv outflows, the SDSS J0838+2955 
result should be divided by 2.5 (since for the Siv outflows we use fl = 0.08). 



6. RELIABILITY OF THE MAIN STEPS IN DETERMINING M and Ek 

Our reported kinetic luminosity (Ek) for SDSS J1106+1939 is an order of magnitude 
larger than the p revious highest value in an established quasar outflow (SDSS J0838+2955, 



Moe et al. 



20091 . see Table [5]). Due to the potential importance of this result to AGN 
feedback processes, we review in this section the steps that were taken in order to arrive at 
this result. Our aim is to address possible caveats or systematic issues that might affect 
this result, especially whether Ek can differ significantly from the most plausible value we 
report in Table E] {Z = 4Zq model). 



6.1. Ionic Column Density Extraction and Their Implications 



Reliable measurements of the absorption ionic column densities (Ajon) in the troughs 
are crucial for determining almost every physical aspect of the outflows: ionization 
equilibrium and abundances, number density, distance, mass flux, and kinetic luminosity. A 
firm lower limit on Nion of a given trough is pr oduced by integratin g the apparent optical 



depth (taod) of the trough across its width (see 



Borguet et al 



2OI2I ). Since outflow troughs 



Table 5. Physical properties of the S iv quasar outflows. 



Object 


log{LBol) 

(ergs s"-*-) 


V 

(km s-i) 


log(C/i/) 


log(7V^) 
(cm-2) 


log(ne) 
(cm~^) 


R 
(kpc) 


M 

(Mo yr-i) 


log(£;it) 
(ergs s"-*-) 


Ek/ LboI 
(%) 


J1106+1139'^ 
J1106+1139'' MF87 
J1106-I-1139 4^0 

J1512+1119'' (C4) 
J1512+1119'' (C2) 


47.2 
47.4 
47.2 

47.6 
47.6 


-8250 
-8250 
-8250 

-1050 
-1850 


o+°-2 

U.U_Q 2 
-U.O_o.2 

-1.7 to -1.4 

-0.9l°:J 


oo Q+0.2 

22.6t°i 

oo 1+0.3 

-^-^•-■■-o.i 
19.8-20.4 

2i.9i2:i 


4 1+0-14 

4 1+0-14 
^■-'—0.37 

^•-■■-0.37 

< 3.3 

5 4+2-7 


O-lSlaoe 
0.29l°i? 

n qO+0.20 

U.O.^_0.14 

>3.1 

0.3-0.01'^ 


11001™ 

iiool™ 

3901?°° 
>3.4 

55-1"^ 


46.41°| 
46.4l°i 
46.01°:? 
>42.1 
43. 8-42. l'^ 


151^° 
101^ 
5lf 

>io-^ 

< 0.02 


J0838+2955 


47.5 


-5000 


-2.01°:^ 


20.8l°J 


3.8 


o 0+1-5 
'J-'-'-l.O 


3001?^° ^ 


45.41°:^ ^ 


0.81°:^ ^ 



a) Solar metalicity model, 
b) The ranges corresponds to the upper and lower error bars on n^, respectively, 
c) Computed using O = 0.2 (see text for discussion). 
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often exhibit non-black saturation, one has to be careful when assessing the actual Nion of 
a given trough. As shown in Paper II, the actual Nion can be 1000 times larger than the 
value inferred from taod- This is the reason why we treat any singlet trough measurements 
(e.g., Sim A1206) as well as heavily blended doublets (e.g.. Civ AA1548, 1551) as lower 
limits. For wider separated doublets (in our case the important P v and Siv), our template 
fitting allows for a clear distinction between two ) The fully saturated case where 

the taod of both doublet components is identical within the measurement errors. In such 
a case only a lower limit can be put on the trough's Nion- b) taod of the blue doublet 
component is significantly larger than the taod of the red doublet component within the 
measurement errors. In that case, an upper limit for the actual Nion is only a few times 
larger than that inferred from T Ann and can be m easured under the partial covering and/or 



power-law absorber models (see 



AravetaL 



20051 ). 



Following these principles, the strength of our 3 most important Nion measurement in 
the SDSS J1106+1939 outflow is as follows: 



• He I*: the two measured troughs show an absorption case that is essentially AOD. 
Therefore, the measurement is very reliable and the small error reflects it's robustness 
(see Table [1]). 

• P V: In our most conservative case, template fitting shows a taod ratio of 1:1.1, which 
due to the quality of the data is still distinguished from the fully saturated case of 1:1 
ratio. The 1:1.1 ratio yields a significantly higher Nion than the AOD case and also 
causes the large associated error bars. 

• S IV: the bottom of the S IV trough is black. Therefore, if our template fitting 
assumption does not hold here, there is a possibility that we severely underestimate 
the true S IV Nion- We note that in such a case the estimated will be larger for two 
compounding reasons: the first being that the Nion ratio of Siv*/Siv will become 
smaller, reducing the deduced ne of the outfiow and hence increasing its distance. As 
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can be see in Equation (3) a larger R yields larger Ek. The second one is that a larger 
total S IV Nion will necessitate a larger Nh in the photoionization solution, which also 
yields a larger value of E^. 



6.2. Photoionization Modelling: Sensitivity to Different SEDs and 

Abundances 



The lack of observational constraints on the incident SED, especially in the critical 
region between 13.6 eV and the soft X-ray region, motivates us to check the sensitivity of 
the abov e results to other AGN SED s. For comparison we use the MF87 SED developed for 



AGN by 



Mathews &: Ferlandl (119871 ). which is both considerably different than the UV-soft 



SED used in our analysis due to it's strong "blue bump" and also used extensively in the 
literature. The best-fit photoionization models are parametrized by log Uh = —0.2 and 
log Nh = 22.6 cm'^ for MF87, within 0.2 dex of the values obtained with the UV-soft SED. 
We note that, as shown in Table 5, the derived mass flow rate and kinetic luminosity using 
the MF87 determined Nh and Uh are consistent with those derived using the UV-soft SED 
to better than 25%. This is due to the functional dependencies of these quantities on Nh 
and Uh, and due to the higher Qh associated with an MF87 SED (see Equations (1) and 
(3)). We conclude that the derived results are only mildly sensitive to other physically 
plausible quasar SEDs. 

In contrast, the derived M and E^ are more sensitive to departure from solar 
metallicity. A comparison of the results in the flrst and third models given in Table 5 shows 
the following: Compared to the Z = Zq model, the Z = 4Zq model has 3 times less M and 



Ek, while at the same time R is 60% larger (due to the lower Uh of the Z 



The photoionizati on reasons for that behavior are discussed in 



Korista et al 



Zf7^ m odel). 



Dunn et al 



(120081 ) and 



(I2OIOI ). We note that since the super-solar metallicity models gives a better 
solution for the outflow, and the upper bound for the metallicity is close to Z = AZq (see 
Section W7L\i . our representative model indeed yields a conservative lower limit on M and 
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6.3. Reliability of the Distance Estimate 

Our distance estimate is derived from measuring tlie ratio of Siv*/Siv column 
densities. If tliis ratio is larger tlian unity tliere is a possibility that the S iv* is actually 
saturated and that the true column density ratio is close to 2. In this case we can only 
establish a lower limit on ng and therefore an upper limit on R. The converse is true for 
cases where the Siv*/Siv column density ratio is smaller than unity, which is the case for 
our SDSS J1106+1939 measurements. Therefore, the distance we derive for this outflow is 
technically a lower limit. We note that using this R value as a measurement instead of a 
lower limit yields a lower limit on Ek (see Equation (2)). Thus, again our reported Ek is 
conservative. This point is somewhat different than the argument we presented in point 3 
of Section 16.11 above as it holds even if the S IV resonance trough is not black. 

6.4. Global Covering Factor (fi) 

As noted in the Introduction, this aspect is one of the major strengths of the current 
analysis. Based on the usual statistical approach for Civ BALs, the range in Q for these 
Siv outflows is between Q=0.2 (the value for all high ionization BALs that are represented 
by Civ troughs) and fi=0.08, which takes into account that only 40% of Civ BALs 
show kinematically-corresponding Siv absorption (see Section [5.2^ . For the energetics 
calculation, we chose the more conservative f2=0.08, which minimizes Ek- 

6.5. Can The Outflows Hide A Signiflcant Amount of Undetected Mass 

The highest ionization species available in our spectrum is O VI for which we can 
obtain only a lower limit for Ajon- This does not allow us to probe the higher ionization 
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material that is known to exist in AGN outflows via X-r ay observations, the so-called Warm 



Gabeletal 



(l2005a] ). the combined 



Absorber (WA) material. For example, from Table 3 in 
column density of the warm absorber in the NGC 3783 outflow is 10-20 times higher than 
that detected in the UV components. At least in one case, UV spectra of outflows from 
luminous quasars show a similar phenomenon where very high ionization lines (NeVIII, 
and Mgx) sho w that the bulk of th e outflowing material is in this very high ionization 



component (see 



Muzahid et al. 



20121 : Arav et al 2012b in preparation). 



It is certainly possible that such a high ionization phase is also present in the SDSS 
J1106+1939 outflow (and in the other two Siv outflows we report here). In that case, the 
total M and Ek could be much larger than reported in Table [5] if the WA material is located 
at the same distance than the UV outflow. 



6.6. Summary: 10^^ ergs s ^ is a conservative lower limit for E^. in the SDSS 

J1106+1939 outflow 

In this section we demonstrated that essentially all possible deviations from our 
assumptions will lead to a higher value of Ek in the SDSS Jl 106+1939 outflow: possible 
higher column density in either Siv or Pv; metallicity lower than Z = AZq] the derived R 
can only be larger (i.e., if the SiV column density is underestimated); we use a conservative 
value for Q, using the one associated with the general high ionization CiV outflows will 
yield E^ larger by a factor of 2.5; the outflow may very well carry a dominant component 
of high ionization material, to which our rest-frame UV observations are not sensitive; in 
that case E^ oc Nh{UV) + Nh(WA), where Nh{UV) is the total hydrogen column we are 
sensitive to using the UV diagnostics in the data analyzed here, and NniWA) is the Nh 
associated with the higher ionization WA material. 
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DISCUSSION 



The powerful BAL outflow observed in SDSS J1106+1939 possesses a kinetic luminosity 
high enough to play a major role in AGN feedback processes, which typically require a 
mechanical energy input of roughly 0.5-5% of t he Eddington luminosity of the quasar 



(IHopkins fc Elvis 



2010 



Scannapieco fc Oh 



200^ respectively). This quasar, being bright 



for its redshift band radiates close to its Eddington limit (i.e. LboI — Ledd)- Therefore, 
with Ek ^ 5%Lboi, it has enough energy to drive the theoretically invoked AGN feedback 
processes. 

How applicable are these results to the majority of quasar outflows? The investigation 
described here gives the first reliable estimates of R, M and for a few high ionization, high 
luminosity quasar outflows. Previously we only had such determinations for low ionization 
outflows, which comprise only ~10% of all observed quasar outflows. Furthermore, the 
absorption spectrum of these two objects looks very similar to the run-of-the-mill BALQSO 
spectra longward of Lya. This is an important point as the vast majority of available 
BALQSO spectra do not extend to wavelengths shorter than llOOA (restframe) and 
therefore do not cover the Siv/SlV* lines. The phenomenological similarity of the CiV 
and Si IV BALs in the objects presented here to the majority of observed BALs suggests 
that a straightforward generalization of the results may be plausible. SDSS J1106+1939 
is the first VLT follow-up observation we obtained from a dedicated search for possible 
Siv/Siv* troughs. In the coming year, we are scheduled to obtain several more X-shooter 
observations of such candidates and will be able to shed more light on this issue. 

The distances found for the three outflows we report here range from ~ 100 parsecs 
to a few kiloparsecs. These are similar to the distances inferred for outflows in which the 
density diagnostic i s obtained from the study of excited troughs of singly ionized species 



as Fell or Sill (e.g. 



Korista et al. 



2008 



Moe et al. 



2009 



Dunn et al. 



2010), but they are 



3-|- orders of magnitu de further away than the assumed acce leration region (0.03-0.1 pc) of 



line driven winds (e.g. 



Murray et al 



1995 



Proga et al. 



20001). This result is consistent with 
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almost all the distances reported for AGN outflows in the literature. However, the current 
research expands the claim to the majority of high ionization outflows. We conclude that 
most AGN outflows are observed very far from their initial assumed acceleration region. 
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